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Charmonia hadro- and photo-production are analyzed in the framework of the color-octet model taking into 
account higher-order effects induced by initial-state radiation of gluons. We argue that color-octet matrix ele- 
ments obtained from Tevatron data might be reconciled with those extracted from HERA data on inelastic J/ip 
photoproduction. Finally we estimate cross sections for events with muons from prompt J/tp's as a background 
for the B physics programme at LHC experiments. 



1. INTRODUCTION 

Over the last years, charmonium physics has 
recovered a leading role in the challenge of under- 
standing both perturbative and non-perturbative 
aspects of the strong interaction dynamics. The 
surplus of J I ip and ip' resonance production found 
at the Fermilab Tevatron Q| triggered an intense 
theoretical activity as conventional wisdom (i.e. 
the color-singlet model) failed dramatically to re- 
produce experimental data. Such a discrepancy, 
by more than an order of magnitude, led to the 
introduction by Braaten and Fleming j|] of a 
new production mechanism, the so-called color- 
octet model (COM), later cast into a well defined 
framework based on an effective non-relativistic 
field theory (NRQCD) [§. However, although 
those non-perturbative parameters introduced in 
the model and needed to bring accordance with 
experimental data can be rigorously defined in 
the effective theory (as vacuum expectation val- 
ues of effective field operators), the possibility of 
obtaining numerical values from lattice calcula- 
tions seems unfortunately still far. This current 
theoretical inability means in practice that such 
long-distance parameters have to be determined 
from fits to experimental data. 
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On the other hand, there is the founded convic- 
tion that due to the distinct energy scales, char- 
monium production cross section at high trans- 
verse momentum factorizes into a short- and a 
long-distance part, the former calculable from 
perturbative QCD, the latter parameterized by 
NRQCD matrix elements expected to be univer- 
sal - i.e. the same for different production mech- 
anisms and energies. However, numerical values 
of some non-perturbative parameters extracted 
from Tevatron Q and HERA || apparently are 
in disagreement each other, casting some doubts 
on the validity of the COM itself, or at least on 
the above-mentioned universality hypothesis. 

This paper is partially devoted to argue that 
the claim on the inconsistency of the NRQCD 
matrix elements obtained from hadro- and photo- 
production may be premature because the effec- 
tive intrinsic transverse momentum of partons 
has been overlooked. 

Indeed, it was well known long time ago 
that particle production in hadron-hadron colli- 
sions requires the introduction, usually as a phe- 
nomenological input, of a primordial k t of partons 
in addition to their expected Fermi motion 
During this Conference, we have heard Gemme 
and Lazzeroni discuss the need of a large intrin- 
sic kt to reproduce experimental data on D and 
B production. In fact, most of the effective trans- 
verse momentum of partons should be attributed 



to multi-gluon emission by the initial state. In 
the literature sometimes this effect is simulated or 
simply parameterized by means of a convolution 
of the hard process cross section with a Gaussian 
kt distribution whose width is an adjustable pa- 
rameter However, strictly speaking the effect 
is modeled more reliably by using a Monte Carlo 
generation of the real complexity of the initial- 
state radiation. As in our previous analysis of 
charmonium hadro-production || we have used 
PYTHIA g managing multi-gluon emission ac- 
cording to an elaborated algorithm . 

2. EFFECT OF THE INITIAL RADIA- 
TION ON THE p t DISTRIBUTION 

Let us point out that, technically, the inclusion 
of multi-gluon emission by itself hardly changes 
the overall cross section for charmonium pro- 
duction obtained from PYTHIA. However, the 
smearing caused by the initial-state radiation is 
responsible for a migration of events towards the 
high-p t tail of the charmonium distribution, thus 
enhancing the differential cross section in the re- 
gion where experimental points actually lie. In 
figure 1 (see Ref. [jll|) this effect is shown, clearly 
pointing out the need for a correction of the re- 
sults from Cho and Leibovich Q if a more realistic 
description of the underlying hadronic complex- 
ity of the interaction (i.e. initial-state radiation) 
is required. 

In taking into account this pt-dependent K 
factor, new (lower) values for the long-distance 
parameters, and ultimately the color-octet ma- 
trix elements, are obtained from charmonium 
hadroproduction. In Table 1 we present our val- 
ues |ji"Tf determined from older Tevatron data 
|l| together with the new ones |lj,[l3| on J/ip 
hadroproduction, both showing a better agree- 
ment H] with HERA results than before 

3. THEORETICAL UNCERTAINTIES 

There is no space to review all the theoretical 
uncertainties related to the extraction of matrix 
elements from charmonium production which on 
the other hand have extensively been reviewed in 
the literature (see for example Let us only 
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Figure 1. Curves obtained from PYTHIA (not 
fit) including the color-octet mechanism for J/ip 
production at the Tevatron using the same pa- 
rameters and parton distribution function as in 
Ref. [4] . The solid line corresponds to initial and 
final state radiation off and the dotted line to ini- 
tial and final state radiation on. 



mention: 

• Ambiguities in the choice of energy scales, 
charm mass, PDF, ... 

• Partial breaking of the heavy-quark spin 
symmetry in charmonia systems, ... 

• Subleading 0(> v 7 ) contributions, higher- 
order a s effects, ... 

• Possible kinematic effects near boundaries 
of phase space 

• Possible effects due to the surrounding 
medium breaking universality 

Below we comment in some detail on the possi- 
ble effects due to the sensitivity to the kinematics 
of soft gluon radiation in the long-distance evolu- 
tion of the cc pair into final charmonium Jig] . 



Table 1 

Color-octet matrix elements (in units of GeV 3 ) from best fits to old and new Tevatron data on prompt 
J/j> production. (M r J/v> =< 0\Oj^ j 1 S )\0 > + r- < 0|O^( 3 P )|0 > /m 2 .) 



NRQCD matrix element: 


<0|O 8 7/v, ( 3 ft)|0> 




Ref. 
Ref. 


11 

13 


] (GRV, r = 3): 
] (r = 3.5): 


(3.4±0.4)xlCT 3 
(2.73±0.45)xlCT 3 


(6.0±1.2)xl0~ 3 
(5.72±1.84)xlCT 3 



3.1. Finite width effects 

Finite width effects can have consequences 
as Ernstrom, Lonnblad and Vanttinen discussed 
some time ago fl6|| . Indeed, even a small en- 
ergy difference (of order Mv 2 ) between the mass 
of final charmonium and the invariant mass of 
the created cc pair should have an effect on the 
(steep) momentum distribution. Nevertheless, 
the point is that the formation of an intermediate 
colored cc state does not preclude the possibility 
of an energy deficit and not necessarily an excess 
w.r.t. charmonium mass, the latter implying the 
emission of soft gluons as usually considered. In 
other words, the cc pair could get energy from the 
surrounding medium. 

This viewpoint is in fact not surprising since in 
the color evaporation model jl7| based on duality 
arguments, the cross section for charmonium pro- 
duction is written as an integral over the s range 
Am 2 and Am 2 D (i.e. the upper limit equal to open 
charm production). However, charmonium pro- 
duction takes place in a "hot" medium as Hoyer 
pointed out in his talk (for more details see Ref. 
18 1), so sub-thresholds effects likely should have 
an influence. In the improved color evaporation 
model this means that the heavy quark pair can 
get energy from the surrounding color field to give 
open charm production beneath s — Am 2 D [|l9| . 

Analogously in the COM, those dipole tran- 
sitions occurring during the evolution into final 
particles should be associated to absorption of 
soft gluons from the neighborhood as well. For 
example a ccpS"^] pair may absorb two non- 
perturbative gluons becoming a 0(v 2 ) Fock com- 
ponent of the J/ip state or may radiate two soft 
gluons thus with nearly unity overlap with the 
J/ip wave function. The conclusion is that on 
the average one should consider both possibili- 



ties, likely decreasing the net impact of the above- 
mentioned effect on the pt distribution. 

4. LHC PREDICTIONS 

As outlined in the previous section, the extrac- 
tion of NRQCD matrix elements from charmo- 
nium production is not yet fully satisfactory due 
to several theoretical uncertainties. Still more 
analysis has to be done to check the universality 
hypothesis allowing comparisons between differ- 
ent production reactions. Polarization of charmo- 
nium from gluon fragmentation at high p t might 
be decisive to verify the validity of the COM [p0| . 

Nevertheless, under the assumption of the va- 
lidity of the COM to explain charmonia pro- 
duction, long-distance matrix elements obtained 
from Tevatron data should parameterize those 
uncertainties to a large extent. Therefore, in ex- 
trapolating to LHC energies, predictions should 
be reliable if correctly performed. In our case, 
this actually means running PYTHIA with the 
same parameters and options (e.g. radiation on) 
as previously employed in the fit to Tevatron ex- 
perimental points. 

In Figure 2 we show our predictions for prompt 
J/ip direct production at the LHC using the 
CTEQ PDF. A caveat is in order, however. In 
changing the choice for the PDF there are no- 
ticeable changes on the production rate by about 
50% in certain regions of p t (see Ref. Jn]] ) . 

4.1. Muon cross sections from prompt J/ip 
production at the LHC 

We have estimated the cross sections for events 
with muons from prompt J/ip's, passing typical 
pseudorapidity and transverse momentum cuts, 
foreseen in the first and second level triggers for B 
physics at the two LHC experiments ATLAS and 



sources for the B physics programme at LHC ex- 
periments, especially regarding the search for CP 
violation via the golden decay mode Bj—>J/ipK s . 
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Figure 2. Our predictions for prompt J/tp direct 
production at the LHC according to the color- 
octet model implemented in PYTHIA for the 
CTEQ PDF [10]. Dotted line: CSM; dashed line: 
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contributions; dot-dash line: 3 S[ 
solid line: all contributions. No- 
tice that asymptotically the 3 s[ 8 ^ contribution 
becomes dominant, likely implying a transverse 
polarization of charmonium at high p t . 



CMS. We shall denote as /i3, /x6 (/j3/j3, /i6/i3) 
those events with at least one (two) muon(s) 
whose p t is (are) larger than 3 GeV or 6 GeV, 
respectively. 

From our modified version of PYTHIA with the 
COM implemented in, we obtain 

ct(/x3) = 622 nb 

<r(fi6) — 52 nb 

cr(/Lt3/Li3) = 30 nb 

<t(/i6/i3) — 9 nb 

In all cases an additional pseudorapidity cut 
1 77 1 < 2.5 for muons was required on generated 
events. 

Those numbers should be useful to evaluate 
an order-of-magnitude estimate of background 
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